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Abstract 

In the framework of the band theory, we consider two tight-binding models of insu- 
lators. The first one, proposed recently by Taraskin et al, is a translationally invariant 
system, built out of two independent non-overlapping bands of single-particle orbitals that 
are coupled by a weak inter-band hybridization. This kind of insulator exhibits unphysical 
properties: we show, in particular, that the one-body density matrix does not depend on 
the width of the gap between the bands. Consequently, there is no derealization effect 
with increasing metallicity. In the second model there are also two bands. However, 
they are not imposed by construction but are created from a band of single-particle or- 
bitals due to the breaking of the translational symmetry by a periodic potential. These 
bands are separated by a gap for all nonzero values of the unique energy parameter of the 
model. We demonstrate that the one-body density matrix has the same structure as in 
the first model. As a result, the large distance asymptotic formulae derived by Taraskin 
et al in dimensions D = 1,2,3, apply as well, but only for very large gap widths. In 
D = 1 and in the diagonal direction of D = 2 cases, we derive a stronger asymptotic 
formula, valid for all gap widths. The both kinds of asymptotic formulae are composed of 
a dimension-dependent power-law factor and a gap-dependent exponentially decaying fac- 
tor. The latter asymptotics implies that the exponential decay rate vanishes linearly with 
the vanishing gap. In non-diagonal directions, we have found numerically that the linear 
scaling is replaced by the square root one. Independently of the direction, the exponential 
decay rate grows logarithmically with the gap width, for sufficiently large gaps. 
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The rapid progress in computational techniques used for calculating properties of solids 
enabled researchers to use a localized real-space approach to describing such properties. It has 
made possible large-scale calculations based on the density functional theory jT], in particular 
calculating the electronic structure of solids by means of O(N) methods EI]. In all these 
methods it is the one-body density matrix (DM) that is the central quantity. Its decay rate 
determines the degree of locality of all relevant quantities and is decisive for the speed of the 
involved algorithms. This is why one observes a growing interest in localization properties of 
DM in recent years. However, the first result concerned with the rate of decay of DM was 
published by W. Kohn as early as in 1959 jl]. It has been established that in a ID model 
of an insulator the DM decays exponentially with distance. The exponential decay of DM 
was reconsidered by Ismail-Beigi and Arias |H] a few years ago (see also the references quoted 
there for developments in the meantime) in a general, model- independent context, for systems 
described by single-particle orbitals in periodic potentials in arbitrary dimensions. Then, He 
and Vanderbilt have discovered the power-law prefactor mutiplying the exponential factor 
in the asymptotic formula for DM. Finally, Taraskin et al [Zj presented their minimal model of 
an insulator and described, analytically and numerically a power-law decay and an exponential 
decay in dimensions D = 1,2,3. 

In this communication we reconsider the model of insulator proposed in [7j. This is a 
translationally invariant system described by two kinds of single-particle orbitals between which 
the electrons can hop. We demonstrate that it is not a kind of insulator considered in [Sj or 
earlier papers on the subject. For instance, for those values of the energy parameters of the 
model for which there exist two bands separated by a gap, there is no relation between the 
decay rate of DM and the size of the gap. We have succeeded in deriving analytically the large 
distance properties of DM of another system, described by one kind of single-particle orbitals 
and a periodic external potential. This system exhibits two bands separated by a gap for all 
nonzero values of the parameter determining the strength of the potential. Specifically, in the 
large distance behavior of DM elements of this system we have determined the power-law factor 
and the exponential factor. We have found also the scaling law of the inverse correlation length 
versus the vanishing gap. 

The model of an insulator proposed in [7j (TDE insulator) is described by the following 
second-quantized Hamiltonian 

In the above expression, *, j represent the lattice sites of a D-dimensional Bravais lattice, 
while stands for a pair of nearest neighbors on this lattice. The operators {at a^} 

create, annihilate, respectively, a spinless fermion in a single-particle orbital whose bare 
energy is e^, where /i differentiates between the two kinds of involved orbitals. The orbitals 
\i,fj) satisfy the orthonormality condition (j, = Sij6 Utfl . The hopping integrals between 

the orbitals \i,fi) and \j, v) are t^, and they are nonzero only if i,j are nearest neighbors. If 
the periodic boundary conditions are imposed, then the system is diagonalized by passing to 
plane- wave orbitals \k,fi), with the wave vector k in the first Brillouin zone. In order to reduce 
the number of energy parameters in the model, we make the bands symmetric. Specifically, we 
equalize the intraband transfer integrals tn = t 2 2 — t, introduce Ae = 61—62, with t > and 
Ae > for definitness, and shift also the zero of the energy scale to (ei + £2) /2. Moreover, we 
express all the energies in the units of the intraband hopping t. Then, the eigenenergies of the 
system, as functions of the wave vector fc, form the upper, A£, and the lower, A^, bands Ajjjr: 

\± = 2S k ±2[(Ae/4) 2 + r*Sl] 1/ \ (2) 
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Figure 1: Curves of constant gap 5 for D = 2. The thick curve corresponds to 5 = 0; it splits 
the t > 0, Ae > quadrant into two disjoint regions of S > and 5 < 0. 

where r = t^/t is the interband transfer integral in the units of t, r > for definitness, 
and Sfc = \ , t ^ exp {ik{i — j)) stands for the structure factor. Clearly, A^ and are 
mutually symmetric about zero and there are two independent energy parameters: r and Ae. 
Let 8 = mirikXk —maxkXZ be the width of the gap (if 5 > 0) or the overlap of the two bands (if 
5 < 0). In Fig.^ we display the curves of constant 5, in the plane (r, Ae), for D = 2 . It is seen 
that if the bare energies of the two kinds of single-particle orbitals are such that Ae > 8, then 
for any hybridization r > the width of the gap 5 > 0. Otherwise, the hybridization has to be 
large enough. Generally, we need Ae > AD to have a gap for arbitrarily small r. Now, consider 
the TDE insulator, i.e. we set 5 > and fill up completely with spinless electrons the lower 
band. Moreover, the underlying Bravais lattice is assumed to be a D-dimensional simple cubic 
lattice, with lattice vectors expressed in the units of the lattice constant (i.e. the components 
ji, I = 1, . . . , D, of lattice vectors j are integer). Then, the zero-temperature matrix elements 
of DM are given by the zero-temperature averages {o^l ^ a i+r,v) which in the thermodynamic 
limit assume the form 



/7T /*7T 
■■•/ dfcexp(ifcr) [l =f AA~ 1 {k 
-7T J —TV 



( a ll a i + r,l) = ( 2 (^r) /••-/ dkex V (tkr)[l T AA~ L (k)\, (3) 
Kifli+r.a) = -(2(2n) D )- 1 f ■ f dfcexp (ikr) S k A-\k), (4) 

J — 7T J — TV 

where A = Ae/Ar and A(k) = (A 2 + S^) 1 ^ 2 . Clearly, the matrix elements of DM depend on 
the two energy parameters of TDE insulator only through their ratio Ae/r, hence in that part 
of (r, Ae)-plane where 5 > 0, they are constant along the rays with the slope 4A, emerging 
from the origin (0, 0). Consequently, the DM does not depend on the width of the gap between 
the upper band and the lower one. In other words, by following a suitable curve in the region 
5 > 0, we can impose any relation between the parameter A and the gap width 5. For instance, 
along the lines of constant r the parameter A is an increasing function of 5, whose minimum 
is A T , and the limit A — > A+, is equivalent to the limit 5 — > + . Only for a sufficiently strong 
hybridization, r > 1, the minimum A T attains zero and then A — > + is equivalent to 5 — > + . 
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The latter limit has been studied numerically in [Zj in order to obtain the large distance behavior 
of DM for 5 — > + . 

Moreover, despite the hybridization of the two types of bands, in some lattice directions there 
are no density-density correlations between the electrons, that is the correlations (tij aj til i) 
vanish identically in some lattice directions. 

Now, consider the model, where only one sort (/x fixed) of single-particle orbitals \j, /i) = \j) 
is present, whose Hamiltonian is 

H cb = Uidfcii + t^2 ( a i a j + h - c ) > ( 5 ) 

where we used the same notation as in (|T} but with the fixed index \x suppressed, and Uj 
denotes a periodic external potential. Suppose that the underlying Bravais lattice consists of 
two interpenetrating sublattices that differ by a primitive translation, i.e. the nearest neighbors 
of a site on one sublattice belong the other one. Then, we set Uj = U\ on one sublattice 
and Uj = U% on the other one. Under periodic boundary conditions the Hamiltonian (jSj) is 
diagonalized by plane wave orbitals \k) with the wave vector in the first Brillouin zone of a 
sublattice. Specifically, shifting the energy scale to (U\ + U-i)j2 and expressing all the energies 
in the units of the transfer integral t, we obtain the upper A^ and the lower bands of 
eigenenergies labeled by k: 

A± = ±2y / ( M /2) 2 + SI = ±2A cb (k), (6) 

where u = (U2 — U\)/2t is the unique energy parameter of the model. Without any loss of 
generality we can set u > 0. As in the TDE insulator, the two bands are mutually symmetric 
about zero and are separated by a nonzero gap 5 cb = 2u, for any nonzero u. For 11 ^ and 
completely filled lower band, the system given by (jSJ) is an insulator in the sense of the band 
theory, which we call the chessboard insulator. In contrast to TDE insulator, the chessboard 
insulator belongs to the class of insulators, for which a definite scaling of the exponential decay 
rate versus the vanishing gap is expected 5J. Let the underlying lattice be a D-dimensional 
simple cubic lattice. For definiteness, we set Uj = U 2 at the sites of the even sublattice. The 
zero temperature, non-diagonal elements of DM are given by 

(ala i+r ) = (2n)- D [ dk exp (ikr) lk 2 {u 2 + {-1)"% - (-1)^ [1 + (-1)°*] «&} , (7) 

J B.Z. 

where the D-dimensional integral is taken over the first Brillouin zone of a sublattice, r ^ 0, 
and the following abbreviations have been used 

D 

a r = J2 r h Pk = 2S k + 2A cb (k) , 7^ = 4p k A cb (k). (8) 
1=1 

Without any loss of generality, we can assume that for all I — 1, . . . , D, ri > 0. Then, a r is a 
(noneuclidean) distance between the lattice point r and the origin. To reveal the structure of 
DM given by (J7J), it is convenient to rewrite it separately for a r even and odd: 

(ata i+r ) = -(2n)- D I dk exp (ikr) S k A^(k), if a r = 2m + 1, (9) 

J B.Z. 

(aja i+r ) = -(-iy i {2(2n) D )- 1 u dfcexp (ikr) A^ ) 1 (fe), if a r = 2m. (10) 

J B.Z. 
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It appears that the large distance behavior of DM elements for the chessboard insulator is 
determined by the function lZ cb (r), 



n cb (r) = (2tt)- d f dfcexp(ifcr)A^(fc), (11) 

J B.Z. 



since, 

1 D 

(afa i+r ) = -- SiK cb (r), if a r = 2m + 1, (12) 
i=i 
u 

(aja i+r ) = --{-lpn cb (r), if a r = 2m, (13) 

where Si1Z c b(r) = Tl c b{j\ + 1, . . . , ro) + TZ c b(ri — 1, ... , m), etc. Note that the matrix elements 
of DM depend on lZ cb (r) evaluated only at the points r with o v even. Similar relations hold in 
the case of TDE insulator. On introducing 

/TT /»7T 
•••/ dfeexp^fcrOA-^Jfe), (14) 
-TT J —TT 

we obtain 

1 D 

(aJ 1 a i+r , 2 ) = --^W(r), (15) 
i=i 

( a tfi + r,l) = TAll(r). (16) 
In what follows, we limit our analysis to the D = 1, 2 cases. Then, one finds easily that 

nr) = {° v<) . i ^- isodd - (17) 
[7c c6 (r)|« = A, if cr r is even . 

Therefore, as in the case of the chessboard insulator, the DM elements depend only on lZ(r) 
evaluated at the points with a r even. Consequently, the matrix elements (a^Oi+r^) vanish 
at the lattice directions with a r even, while (a^ ^cii +rai ) vanish at the lattice directions with 
a r odd. Moreover, the large distance behavior of the DM elements of both insulators under 
considerations, is given by the function TZ c b(r) restricted to the points with a r even. In all the 
cases studied here we have found that up to a coefficient independent of a r 

K cb (r)~a^eM-°r/t;(e)), (18) 

for sufficiently large cr r , where e stands for an energy parameter and £(e) is the correlation 
length. 

Apparently, a straightforward way to obtain a large a r asymptotics of lZ c b{r) is to expand 
A cb (k) in powers of u~ 2 , for large u, carry out the integrals of products of cosine functions 
and to approximate the Euler T-functions, that arise, by the leading term of the Stirling's 
asymptotic expansion. Then, in D — 1 case we obtain the inverse correlation length 

QeH = ln« + «" 2 , (19) 
while in directions of finite nonzero slope x = r i/ r 2 of the D = 2 case, 

C» = ln« + u~\2 + X + x- 1 ) ~ £ (1 + XT 1 ln(l + X°), (20) 

Q=-l,l 
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Figure 2: The calculated numerically exactly function lZ c b(r), denoted p nume r, divided by its 
TDE asymptotics, denoted p aS ympt, which up to a constant is given by (|2*0jl. versus distance a r , 
for D — 2, % = 1, 2, 4, and two values of it. 

and the exponent 7 = D/2. A similar procedure gives asymptotic formulae along the axes. 
The asymptotics (jTSjl with fTH|) or l}2Uj) - the TDE asymptotics, on setting w = 2A amounts to 
the corresponding asymptotics derived in [Jj. 

While the TDE asymptotics obtained for the chessboard insulator and the coresponding 
one for the TDE insulator have the same structure, their physical content is very different. 
The energy parameter u is in one to one correspondence with the gap width of the chessboard 
insulator. In contrast, A is in no definite relation with the gap width of the TDE insulator. 

To test the quality of the asymptotic formula (j!8)) with (|2()jl. one can compare it with a 
calculated numerically exactly function lZ c b(r). The log-log plot of both functions is to rough 
for such a test. We have plotted, in Fig. El TZ c b(r) calculated numerically exactly, divided 
by its TDE asymptotics. It is seen that, in the case of D=2 chessboard insulator, the TDE 
asymptotics is reasonably good only for as large values of u as 16 or higher. 

To obtain an asymptotics stronger than the TDE asymptotics, valid for small and large 
values of u, we have applied to the integrals in (jllj) the Laplace asymptotic expansion. Unfor- 
tunately, we have succeeded to carry out the calculations only in the D = 1 case (more results 
on the correlation function in this case can be found in |8 ) and for the diagonal direction in 
the D = 2 case. The result is of the form (|18j) . with the inverse correlation length 

r X (e) =He + V^Tl), e=(2D)~ 1 u, (21) 

and the exponent 7 = D/2. Making a similar comparison as that displayed in Fig. El we find 
a satisfactory agreement, see Fig. El Apparently, in the limit of the vanishing gap the inverse 
correlation length (|21j) vanishes proportionally to the gap width. In the one-dimensional case, 
the same result has been obtained numerically for periodic potentials of period higher than 2 
Q. 

Assuming that the asymptotic behavior of TZ c b(r) is of the form (|T5|) with the exponent 7 
independent of u, we have found numerically, see Fig. |H, that in non-diagonal directions the 
inverse correlation length vanishes as y/u as u — > 0. For sufficiently large u, depending on the 
direction, the correlation length behaves as in the diagonal case, i.e. obeys (|2*T|) . Our data show 
also that the decay rate of DM elements is the slowest in the diagonal direction. 
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Figure 3: The calculated numerically exactly function lZ c b(r), denoted p„ umer , divided by its 
asymptotics p asy mpt that up to a constant is given by (}2*T]) . versus distance a r , for D = 2, x — 1 
and three values of u. 




Figure 4: The inverse correlation length, obtained from the calculated numerically exactly 
function TZ c b(r), versus u, for D = 2 and three non-diagonal directions. Thin continuous curves 
are the fits by the function consty/u. The thick curve represents £ _1 in the diagonal direction, 
according to (J5TJ). 
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In summary, this work has been inspired by the recent paper of Taraskin et al [7] on spatial 
decay of DM in a translationally invariant model of insulator - the TDE insulator. Having 
found that the TDE model shows unphysical features, we have carried out similar studies of 
another model of insulator - the chessboard insulator, whose translational symmetry is broken. 
In contrast to TDE insulator, in the chessboard insulator there is a unique energy parameter 
which determines both, the gap between the bands and the DM elements. We have succeeded 
in deriving analytically a dimension-dependent power-law factor, predicted in j^J, with the 
exponent 7 = D/2, and a gap-dependent exponential decay rate, discussed extensively in 
jHj- Specifically, in D = 1,2,3 cases, the DM elements of the chessboard insulator obey the 
asymptotics derived for TDE insulator by Taraskin et al, which we have found to hold only 
for very large gap widths. In D = 1 case and along the diagonal in D = 2 case, we have 
derived a stronger asymptotic formula, valid for arbitrary gap widths. The latter asymptotics 
implies that the inverse correlation length of the chessboard insulator vanishes linearly with 
the vanishing gap width. Concerning non-diagonal directions, on assuming that the exponent 
7 is independent of the gap width, we have found numerically that the inverse correlation 
length vanishes as the square root of the gap width. Independently of the direction, the inverse 
correlation length grows logarithmically with the gap width, for sufficiently large gaps. A more 
comprehensive discussion will be given in a separate article jH]- 

T.K. is grateful to the Institute of Theoretical Physics of the University of Wroclaw for kind 
hospitality and financial support. 
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